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Abstract 

We consider a model of n non-intersecting squared Bessel processes with one 
starting point a > at time t = and one ending point 5 > at time t — T . After 
proper scaling, the paths fill out a region in the te-plane. Depending on the value of 
the product ab the region may come to the hard edge at 0, or not. We formulate a 
vector equilibrium problem for this model, which is defined for three measures, with 
upper constraints on the first and third measures and an external field on the second 
measure. It is shown that the limiting mean distribution of the paths at time t is 
given by the second component of the vector that minimizes this vector equilibrium 
problem. The proof is based on a steepest descent analysis for a 4 x 4 matrix valued 
Riemann-Hilbert problem which characterizes the correlation kernel of the paths at 
time t. We also discuss the precise locations of the phase transitions. 



1 Introduction 

Let X{t) be a squared Bessel process with parameter a > —1, i.e., it is a diffusion 
process (a strong Markov process with continuous sample paths) on R"*" with transition 
probability density given by (cf. [35] ) 
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is the modified Bessel function of the first kind of order a; see [1^, p. 375]. If d = 2(a + 1) 
is an integer, the squared Bessel process can be viewed as the square of the distance 
to the origin of a d-dimensional Brownian motion. Some apphcations of this diffusion 
process can be found, for example, in [29j. 

We consider a model of n independent copies Xj{t), j = 1, . . . ,n, conditioned such 
that for some T > 0, 

• Xj (0) = aj,Xj (T) = bj for J = 1, . . . , n for some given values < oi < a2 < • • • < 
On, < 6i < 62 < • • • < 6„, 

• the paths do not intersect each other for < t < T. 

It follows from the seminal paper of Karlin and McGregor [50] that the positions of 
the non-intersecting squared Bessel paths at any given time t G (0, T) are a determinantal 
point process. This means that there exists a correlation kernel Kn such that for each 
positive integer m the m-point correlation function takes the determinantal form 

(1.4) 

Indeed, we have by |30) that the positions of the paths at time t have a joint proba- 
bility density 

VnA^l, - ■ ■ ,Xn) = -;^<iet[fj{xk)]j,k=l,...,ndet[gj{Xk)]j,k=l,...,n, (1-5) 

where fj{x) = pf{aj,x), gj{x) = p^_^{x,bj), j = 1, . . . ,n, and is the normalization 
constant 



Zn,t= Vn,t{xi,...,Xn)dxi...dXn. (1.6) 

^{0,00)" 

Hence, (jl.Sp is a biorthogonal ensemble [8], which is a special case of a determinantal 
point process. The correlation kernel is given by 

n 

Kn{x,y) = ^(t)j{x)i}j{y), (1.7) 
i=i 

where the functions (pj, ^j, j = 1, . . . ,n are such that 

span((/)i, ...,(/)„)= span(/i, . . .,/„), span(V'i, ...,ipn)= span(5ti, . . . ,5„), 

and 

(j)i{x)ipf,{x)dx = 5i^k, i,k = l,...,n. (1.8) 







We are interested in the situation where all paths start at the same position at t = 
and end at the same position ai t = T. Different types of initial and ending conditions 
are considered in [3H [33| I34j . among others. In the confluent limit aj — )■ a > and 
6j — 7- 6 > 0, the biorthogonal ensemble structure is preserved. 
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Figure 1: 50 rescaled non- intersecting squared Bessel paths with a = b = 2 (left) and 
a = 6 = 1/4 (right). 



If a = 6 = 0, i.e., all paths start and end at the origin, this system is known as chiral 
or Laguerre ensemble in random matrix theory, cf. \19\ [32} [35l 147) . The case where 
a > and 6 = was considered in [S^ and In that situation, all paths, after proper 
scaling, initially stay away from the hard edge at x = 0, but at a certain critical time t* 
the lowest paths hit the hard edge and are stuck to it from then on. The limiting mean 
density of the paths at time t is characterized by a non-real solution of an algebraic 
equation of order three; see also [37| Appendix] and 31] for an interpretation in terms 
of a vector equilibrium problem with two measures. 11 1 ^ t* , the local correlations obey 
the usual scaling limit from the random matrix theory, leading to the sine. Airy, and 
Bessel kernels [37]. A new kernel was established near the origin at the critical time t* 
in [38]. The kernel admits a double integral representation which resembles the Pearcey 
kernel [ZlEllII]. 

In this paper, we consider the case where both a and b are positive. To obtain 
interesting results, we make a time scaling as in [37J so that the time variable depends 
on the number of paths n. That is, we replace the variables t and T by 

so that < t < 1. Now, letting n — )• oo, the paths fill out a region in the tj;-plane that 
looks like one of the regions shown in Figures [1] and [21 depending on the product ab. 
Suppose a and b stay far away from the hard edge x = 0, such that ab > 1/4. Then we 
are in a situation as illustrated by the left picture in Figure [TJ All paths remain positive 
for all time t. On the other hand, if ab < 1/4, there are two critical times ti and t2, 
as depicted in the right picture of Figure [TJ All paths initially stay away from the hard 
edge (a wall) and the lowest paths hit the hard edge at the time ti. Since all paths are 
required to end at a positive position b, the lowest paths are stuck to the hard edge only 
for a while, and they leave the wall at a second critical time t2. The times ti and t2 will 
come together if ab = 1/4, and we are then led to an intermediate situation shown in 
the left picture in Figure [2j The paths fill out a region bounded by the solid line, which 
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Figure 2: 50 rescaled non-intersecting squared Bessel paths with a = 6 = 1/2 and 
= 1/2 (left) and around the point tc (right). 



is tangent to the hard edge at a multicritical time tc', see also the right picture in Figure 
[2] for a magnified view around tc- 

It is the aim of this paper to describe the region that is filled in the large n limit 
by the non-intersecting squared Bessel paths with one positive starting point and one 
positive ending point. Our results are stated in the next section. 



2 Statement of results 

2.1 A vector equilibrium problem 

The main ingredient of our analysis is a vector equilibrium problem involving both 
an external field and two upper constraints. To this end, we define, as usual (see |46j), 
the logarithmic energy of a measure fi by 

I{fi) = // log ^^d//(x)d/i(y), (2.1) 
JJ \x-y\ 

and the mutual energy of two measures /i, v by 

/(/i, i^)= [[ log -^dfi{x)diy{y). (2.2) 
J J \x-y\ 

Consider the energy functional 

3 2 „ 

E{ui,U2,i^3) = ^I{i^j)-Yj^^'^j^'^j+^^+ y{x)djy2{x), (2.3) 
i=i i=i 

where 

Vix) = _2Vax_2Vb^ 
^ ^ t{l-t) t l-t ^ ^ 

The vector equilibrium problem is to minimize E{vi,V2,i'-i) among all measures i^i, 
1^2 and z^3 satisfying the following conditions (E1)-(E3): 
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(El) vi is a measure on M := (— oo,0] with total mass 1/2 that satisfies the constraint 

ui < pi, (2.5) 
where pi is the measure on M~ with density 

|x|~^/^ xGM". (2.6) 



dpi _ Va|^|„i/2 



dx irt 



(E2) 1^2 is a measure on := [0, oo) with total mass 1. 

(E3) f3 is a measure on M~ with total mass 1/2 that satisfies the constraint 

1^3 < P3, (2.7) 
where p^ is the measure on with density 

dp3 Vb ,^^,_i/2 



dx 7r(l — t) 



x\'^'\ xeR'. (2.. 



An electrostatic interpretation of the above equilibrium problem is the following: 
Consider three types of charged particles. The particles of the first and the third types 
are put on R~, both having total charge ^. The particles of the second type are put on 
R"*" with total charge 1. Particles of the same type repel each other. The particles of 
the first and the second types interact with each other with an attraction that is only 
half of the strength of the repulsion of particles of the same type. So do the particles 
of the second and the third types. The particles of the first and the third types do not 
interact directly. The particles of the second type are influenced by the external field V. 
The measures pi and ps in ()2.6p and ()2.8p play the role of upper constraints for the first 
and the third measures. Namely, we minimize the equilibrium problem with the extra 
condition that the particle densities of the first and the third types do not exceed the 
densities of pi and ps, respectively. 

Vector equilibrium problems with mutual interaction can be traced back to \26\ [27] , 
see also [l5]. The equilibrium problem with constraints appeared before mainly in the 
asymptotic analysis of discrete orthogonal polynomials, see e.g. [5} f2D]. 

Our first result concerns the structure of the minimizer of the vector equilibrium 
problem. 

Theorem 2.1. Let a > 0, 6 > and < t < 1. Then there is a unique minimizer 
{pi, p2, fJ-s) of the energy functional ()2.3p under the conditions (El), (E2), (E3). The 
minimizing measures satisfy: 

(a) The support of pi is the negative real axis and there exists ri > such that 

supp(pi - pi) = (-00, -ri]. (2.9) 
The density of the measure pi — pi vanishes like a square root at — ri if ri > 0. 
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(6) The support of fi2 is bounded and there exist < p < q such that 

supp(^2) = [p,g], (2.10) 

and ji2 has a density which is real analytic in the interior of the support and vanishes 
like a square root at the endpoint q and at p if p > 0. 

(c) The support of fi^ is the negative real axis and there exists r^ >0 such that 

supp(/)3 - /Us) = (-cx),-r3]. (2.11) 

The density of the measure — /U3 vanishes as a square root at —r^ if r-^ > 0. 

Note that p, q, ri, and r^ depend on a, b and t. 

The proof of Theorem 12.11 wUl be given in Section O Throughout the whole paper, 
we win denote 

Ai := (-00, -ri), A2:= {p,q), and A3 := (-00, -ra), (2.12) 

i.e., Aj, j = 1, 2, 3 is the interior of the sets supp(/Oi — /^i), supp(^2) and supp(/93 — fi^), 
respectively. 

2.2 Classification into cases 

It turns out that at least one of ri, p, r^ is equal to 0, and generically no two 
consecutive ones are 0. Generically we also have a 0(|a;|~^/-^) behavior of the density at 
the hard edge 0. Then there are three generic cases. 

Case I ri = r3 = and p > 0. The densities of pi — pi and p^ — p^ blow up like 
0(|x|~^/2) as X ^ 0-. 

Case II p > and one of ri,r3 is zero and the other is positive. Thus there are two 
subcases which play a symmetric role. 

subcase Ila ri > 0, p > and r^ = 0. 
subcase lib ri = 0, p > and r^ > 0. 

Case III ri > 0, p = and r^ > 0. 

In Case I we have to add the condition that the densities of pi — pi and /93 — p^ 
blow up like an inverse square root at in order to be in a generic case. There is also a 
non-generic case where one of these densities vanishes as a square root. This is exactly 
what happens in the transition from Case I to Case II. See Case V below. 

The non-generic cases are: 

Case IV Exactly two consecutive values from ri,p, r2 are zero. Thus p = and one of 
ri,r3 is also zero, while the other is positive. 
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Case V ri = r3 = and p > 0. One of the densities of pj — fij vanishes hke a square 
root 0(|x|i/2) 

— )• 0—, while the other blows up like 0{\x\ ^/^) as x 0—. 
Case VI All three are zero: ri = p = r-^ = 0. 

The Case I corresponds to the situation as shown by the left picture in Figure [TJ In 
this case ab > 1/4 and all paths stay away from the hard edge at x = for all t £ (0, 1). 
Both constraints are never active. 

The right picture of Figure [T] which deals with ab < 1/4, contains the three Cases II, 
III, IV. If t < ti OT t > t2, then we are in Case II, where only one constraint is active. If 
ti < t < t2, then we are then in Case III, where both constraints are active. The times 
ti and t2 correspond to Case IV, which represents a phase transition between soft edge 
and hard edge. 

The final Cases V and VI are contained in Figure [2] for the situation ab = 1/4 when 
the domain just touches the wall. We are in Case V if t 7^ tc, and in case VI if t = tc- 

The Cases IV-VI are non-generic critical cases. They correspond to certain phase 
transitions between the generic cases; see also Section 14.41 for a phase diagram in the 
situation where ab < 1/4. 

We shall only consider the generic cases in this paper. For convenience (and without 
loss of generality), it is assumed that 

rs > ri (2.13) 
in Cases II and III, which implies that A3 C Ai C holds in all three cases. 



2.3 Variational conditions 

Denote by 

U^{x) = [ log-^dfi{y) (2.14) 
J \x-y\ 

the logarithmic potential of a measure fi. The minimizer {pi, fJ's) is characterized by 
the following Euler-Lagrange variational conditions: 

Proposition 2.2. The measures pi, ^2 cLnd fi^ satisfy for some constant I G M, 

2U^'^{x) = U^''{x) + Uf'^-{x) -V{x) + l, xEsupp(;U2), (2.15) 
2U^'^x)>U^''{x) + U^'''{x)-V{x) + l, X G M+\supp(/i2), (2.16) 

where V{x) is given in l\2A^ : 

2[7'^i(x) = U^'^ix), X G supp(/9i - ^1), (2.17) 

lU^'^ix) KU^'^ix), X G C\supp(/3i -/ii), (2.18) 

and 

2U^''{x) = U^'^{x), xGsupp(/^3-^3), (2.19) 

2U^'^x) kU^'^x), xGC\supp(/93-/U3). (2.20) 
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Note that the inequahties in (j2.16|) . (|2.18|) and (|2.2U|) are all strict. The strictness of 
the inequality in ()2.16p follows from the strict convexity of the external fields acting the 
second measure; see Lemma 13.31 below. 

The proof of Proposition 12.21 will be given in Section [3l 

2.4 Limiting mean distribution 

The following is the main theorem of this paper. 

Theorem 2.3. Let a,b > 0, < t < 1, and let Kn denote the correlation kernel for 
the positions of the non-intersecting squared Bessel paths at time t in the resettled time 
variables ()1.9p . Let (^1,^2,^3) be the minimizer of the vector equilibrium problem stated 
in Section \2.1\ Suppose that we are in one of the generic cases I, II, or III, described 
above, then 

lim -KJx,x) = ^(x), X e M+. (2.21) 
n->oo n dx 

The limit ()2.2ip expresses that the density of the measure /i2 is the limiting mean 
density of the non- intersecting squared Bessel paths at time t as n — )• 00. 

It is an unusual feature that it is the second measure which plays the main role here. 
In recent papers [22\ [23] concerning eigenvalue distribution of a two-matrix model in 
random matrix theory, a vector equilibrium problem with three measures appears as 
well. In that situation, however, the first component of the minimizer is relevant to the 
limiting mean distribution. 

We prove Theorem (23] by means of a steepest descent analysis of a relevant Riemann- 
Hilbert problem. From this analysis, we also obtain the universality results for the local 
path correlations. For Cases I, II, III, we obtain the sine kernel (in the bulk). Airy 
kernel (at the soft edges) and Bessel kernel of order a (at the hard edge 0), which are 
the typical limiting kernels in random matrix theory, see e.g. |14[ [25t I43| . We will not 
discuss this any further but instead refer to the papers [H El [3 [22] for a more detailed 
analysis in similar situations. 

We expect the conclusion of Theorem 12.31 to remain valid in the Cases IV- VI. The 
local path correlations at the origin in Case IV should be the same as those derived in 
[38] , since they both correspond to phase transitions between a soft and a hard edge. The 
Cases V-VI represent new critical phenomena. It is conjectured that they are related to 
the inhomogeneous Painleve II equation, which bears some connections with the 'critical 
separation' case for non-intersecting Brownian motion with two starting points and two 
ending points [21 [T71 [TSl [21] • These results will be reported in a future publication. 

Theorem 12.31 will be proved in Section [6] 

2.5 About the proof of Theorem 12.31 

The proof of Theorem 12.31 relies on the Riemann-Hilbert (RH) problem for multiple 
orthogonal polynomials of mixed type and its connection to the correlation kernel (jl.7p . 
To see this, we start with a proposition which gives us the explicit structure of the 
biorthogonal ensemble in the confluent case. 



8 



Proposition 2.4. In confluent case — )• a > 0, 6j — >• 6 > 0, j = 1, . . . , n, the positions 
of the non-intersecting squared Bessel paths at time t E (0, T) are a biorthogonal ensemble 
(jl.Sp with functions 

f2j-i{x)=x^-^pf{a,x), j = l,...,ni:=\n/2], (2.22) 

f2j{x) = x^^'^p^+^{a,x), j = 1, . . . ,712 := n - ni, (2.23) 

g2j-i{x)=x^-^p^_tix,b), j = l,...,ni, (2.24) 

g2j{x) = x^'^p^Z]ix,b), j = l,...,n2. (2.25) 

Proof. The functions ()2.22p and (|2.23p follow as in the proof of Proposition 2.1 in |37j . 

In the confluent limit bj — )• b, the linear space spanned by the functions x i— >■ 
Px-ti^j j = Ij • • • ) ^) tends to the linear space spanned by 

x^-^—^p^_t{x,b), j = l,...,n. (2.26) 

Using the differential relations satisfied by the modified Bessel functions of the first kind 
(see [1]), we can prove that 

d 1 

Q^P?{x,y) = —{pt~\x,y) -pf{x,y)), 
y^prHx,y) = YtPnx,y) - (I - {a-l))pr\x,y). 

Now it is easy to show inductively that the linear span of (|2.26p is the same as the linear 
space spanned by the functions ()2.24p and ()2.25p . □ 

Next, we introduce two sets of weight functions wi^i,wi^2 and 1,^2,2, 

'2n^/ax 



Wl 



t 



W2,i{x) = X "/^e , W2,2{x) = X ^" ^>''e ^-^ I^-i ^ _^ • 

which are collected in two row vectors 

^1 = {wi,l,Wi,2), = {W2,1,W2,2)- (2.27) 

Note that by ([□]) and (US} 

/2j_i(x) oc x-^'~^u;i,i(x), f2j{x) OC X^'^Wl^2{x), 
g2j-i{x) oc x-^'"^u;2,i(x), g2j{x) oc X^~^W2,2{x), 



(2.28) 



where fj, gj, j = 1, . . . ,n are given in ()2.22p - ()2.25p . The proportionality constants are 
irrelevant for us. 

We then look for a 4 x 4 matrix valued function Y satisfying the following RH 
problem: 
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(1) y is defined and analytic in C \ . 

(2) For X > 0, Y possesses continuous boundary values Y-^{x) (from the upper half 
plane) and Y^{x) (from the lower half plane), which satisfy 



Y+{x) = Y^{x) 



h W{x) 



h 

where 1^ denotes the k x k identity matrix and W{x) is the rank-one matrix 
W{x) = wf (x)w2(x) 

-J Vx^a[ — t — j^a-^y^rrt 
VI 

V J 



(2.29) 



t 



'a \ t 

2ny/ax 



(3) As z — >• oo, z £ C \ M"'', we have 

Y{z) = (/4 + o(i/z)) diag(z"^z"^z-"^z""2). 

(4) Y{z) has the following behavior near the origin: 

Y{z) = O 
as z 0, z G C \ M+, where 

and the O condition in (|2.32p is taken to be entrywise. 



(2.30) 



(2.31) 



A 


1 


h{z) 


h{z)^ 




Ihz) 


Hz) 


Hz) 


h{z 


)\ 


1 


1 


h{z) 


h{z) 


, Y{z)-'=0 


1 


1 


1 


1 




1 


1 


Hz) 


Hz) 


1 


1 


1 


1 




V 


1 


Hz) 


Hz)) 




I 1 


1 


1 


1 


1 



(2.32) 



\zr. 

log \ z\ 



if - 1 < a < 0, 
if Q = 0, 
if a > 0, 



(2.33) 



There exists a unique solution Y to the above RH problem which involves multiple 
orthogonal polynomials of mixed type associated with the modified Bessel weights wi(x), 
W2(a:;) and their Cauchy transforms; see [12] for details. Due to the singularity of the 
weight matrix near the origin, we need the condition (|2.32p to ensure the uniqueness of 
the solution. 
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It also follows from the results in [12J that the correlation kernel (|1.7|) admits the 
following representation in terms of the solution of the RH problem: 



Kn{x,y) 



1 



2-Ki{x — y) 



(o W2,i{y) W2,2{y))Y+^{y)y+{x) 



Wl,2(x) 



V 











(2.34) 



The representation (|2.34|) is based on the Christoffel-Darboux formula for multiple or- 
thogonal polynomials of mixed type, see also [3]. 

We will obtain the limit (|2.21|) from a Deift/Zhou steepest descent analysis [T ^ ITH t fTE] 
for the RH problem for Y . It consists of a series of explicit and invertible transformations 



(2.35) 



which leads to a RH problem for a matrix valued function K which is such that K tends 
to the identity matrix as n — )• oo. Analyzing the effect of the transformations ()2.35p we 
obtain the hmit (g^I]). 

The rest of this paper is organized as follows. In Section [3] we analyze the equilibrium 
problem and establish Theorem 12.11 and Proposition 12. 2i In Section [5] we introduce the 
Riemann surface built from the solution to the vector equilibrium problem and define 
some auxiliary functions. Section [5] contains the steepest descent analysis of the RH 
problem for Y{^z) and Section [6] gives the proof of the main Theorem 12.31 Throughout 
this paper, it is also assumed that n is an even number so that 

m = n2 = n/2. (2.36) 

This assumption is not essential and is only made to simplify the presentation. 



3 Analysis of the vector equilibrium problem 

The main purpose of this section is to prove Theorem 1 2 . 1 1 and Proposition 12.21 which 
yields the existence and main properties of the minimizers ii\ , ^2 and /i3 . Our basic tool 
is the balayage (sweeping out) of a measure. 

Let be a finite positive measure on a closed set K with positive capacity. We 
denote with \\iy\\ the total mass of z^. A positive measure z) on is called the balayage 
of ly if ||i^|| = and 

J7^(x) = U^{x) + C, q.e. x £ K, (3.1) 

for some constant C, where "q.e." means quasi-every, i.e., with the exception of a set 
with zero capacity, and is the logarithmic potential (|2.14p . see [461 for background 
on potential theory. To emphasize the dependence on K, we also write 

i' = Bal{u,K). (3.2) 
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A relation between v and Bal(z^, K') is established as follows: 

Bal{i^,K) = j Ba\{6s,K)di^{s), (3.3) 

where 5z denotes the Dirac delta measure at the point z. 

For the purpose of this paper, we are interested in the case where 

K = Kc = (—00, — c], with c > 0. 

Then C = in (j3.ip . since Kc is unbounded. By using simple contour integration, it 
can be shown that, if s € the balayage of 5s onto Kc has the density 



~A = T, X e ^c- (3.4) 

dx vr y/|x + c|(s - x) 

In particular, for c = 0, 

dBal(,5„M-) ^. dBal(,5„Ke) 1 ^/i ^ ^_ 

— lim = = , X G M . (3.5) 



dx c^o+ dx vr yT^(s - x) 

Observe that the density (j3.4p is decreasing as |x| increases, with a decay rate 0(|x|~^/^) 
as |x| — )• 00. It will become clear in what follows how the concept of balayage measure 
is linked to the equilibrium problem. 

3.1 Equilibrium problems for v\ and 1^3 

If V2 and are given, the equilibrium problem for v\ is to minimize 

/(z^) lJ''^{x)dv(x), (3.6) 

among all measures v on M~ with total mass 1/2 and satisfying the constraint v < 
where p\ is the measure on M~ with density given by (|2.6p . Clearly, this problem only 
depends on V2. The properties of v\ are given in the following lemma. 

Lemma 3.1. (The measure v\) Suppose 1^2 oind are fixed so that conditions {E2) 
and {E3) of the equilibrium problem are satisfied. Then the measure vi that minimizes 
E{h'i,i^2,'^3) subject to the condition (El) exists. It has the following properties: 

(i) If the constraint pi is not active, then the measure vi is given by 

Z.1 = ^Bal(z.2,K-), (3.7) 



or, equivalently, in view of ()3.3p and ()3.5p . 

di'i{x) 1 f ^/s 



dx 2TXxf\x\ 



dv2{s), xGM". (3.8) 
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(a) If the constraint pi is active, then there exists a constant ri > 0, such that 

supp(/9i - i^i) = (-00, -ri]. (3.9) 

Moreover, we have 

hu''^{x) = U''^{x), X G supp(pi - z^i), 
[lU^^ix) kU^^x), X G R- \supp(/9i - i/i), 

and pi — vi vanishes like a square root at —ri. 

{Hi) Suppose U2 has a compact support. Then here exists a constant K > depending 
only on maxsupp(f2) such that 

Proof. (i) Suppose that vi = Bal(i/2, IR")/2. Then 2f/''i(x) = V^^ix) for x G M", 
which is the variational condition for ()3.6p without constraint. This proves part 

(i)- 

(ii) Part (ii) is proved by using iterated balayage algorithm of [36] . The main idea is to 
construct a sequence of real numbers {ri ^(t)} and a sequence of measures {z^i,fc} 
with limit ri(i) and 1^1, respectively. Since this argument is similar to those given 
in [22J, we omit the proof here. 

(iii) Suppose 1^2 has a compact support and let 

b = maxsupp(f2). 

By part (ii), the constraint is active on some interval [— ri,0) (which could be 
empty). Since ||z^i|| = 1/2, we then have that ri < ci where ci > is the value 
for which pi{[—ci,0]) = 1/2. By properties of balayage we then have that the 
restriction of ui to Kc^ is less than the balayage of 1^2 /"^ onto Kcj^. Hence by (j3.3p 
and (13. 



duiix) 1 f a/F+ci , , , , -, , , 

— ^ < ,, / -di^2{s), X G -00, -ci . 3.12 

dx 27ry^\x + ci\ J s-x 

The function s 1— >• is increasing for s G [— ci, —x — 2ci]. Thus if x < — 2ci — b, 

then the maximum on supp(z^2) is taken at its right end point s = b, and therefore 



— ; < , — , X G (— 00, — 2ci — . (3.13) 

dx - 2^v^|x + ci| b-x ' V , i J V ; 

Then there is a constant K, which only depends on b, such that ()3.1ip holds for 
X < — 2ci — b. Since z/i < pi, we can then adjust the constant if necessary, such 
that ()3.1ip holds for every x < 0. 

□ 
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If z^i and 1^2 are fixed, tlie equilibrium problem for is to minimize (|3.6|) among all 
measures on M~ with total mass 1/2 and satisfying the constraint i' < ps, where 
is given in (j2.8p . Hence, the properties of 1/3 can be derived in a manner similar to the 
analysis given above. We omit the proof and state the results in the following lemma: 

Lemma 3.2. (The measure 1^3) Suppose vi and o-'^^ fixed so that conditions {El) 
and {E2) of the equilibrium problem are satisfied. Then the measure that minimizes 
E{i'i,i'2,'^3) subject to the condition {E3) exists. It has the following properties: 

(i) If the constraint ps is not active, the measure is given by 

[ii) If the constraint is active, there exists a constant r^lt) > 0, such that 

supp(/93 - 1/3) = (-00, -rsit)]. (3.15) 

Moreover, we have 



2U''^{x) = U^^{x), XGSUpp(/93-'^3), 
2C/^3(x) < V'^ix), X G R- \ supp(/93 - 1^3), 



(3.16) 



and p3 — z^3 vanishes like a square root at — r3(t). 



{Hi) Suppose V2 has a compact support. Then there exists a constant K > 0, depending 
only on maxsupp(z^2) such that 

dv-^ix) K ^ , , 

3.2 Equilibrium problem for z/2 

If vi and 1/3 satisfying conditions {El) and (£'3) are fixed, the equilibrium problem 
for V2 is to minimize 

I{v) + j {V{x)-U'''{x)-U'"'{x))dv{x) (3.18) 

among all probability measures v on M^, where V is given by ()2.4p . This is a usual 
equilibrium problem with external field 

V{x)-U''^{x)-U''^{x) (3.19) 

see e.g. [li^ 
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Lemma 3.3. (The measure U2) Suppose vi and are fixed so that conditions (El) and 
{E3) in the equilibrium problem are satisfied, in particular we have Vj < pj for j = 1,3. 
Then the measure z/2 that minimizes i?(zvi, z/2, z^s) subject to condition {E2) exists and 
has the following properties: 

(i) The external field ()3.19p is strictly convex on M"*", and there exist two constants 
< p < q such that 

supp{u2) = [p,q]. (3.20) 

(ii) If ^ supp(pi — h'l) U supp(/93 — z/3) (so that both constraints are active in a 
neighborhood ofO), then ()3.19p is real analytic on M"^ (including 0). 

Proof. (i) The second derivative of (|3.19p is, for x > 0, 

V"{x) - (C/^i)"(2;) - {U'^'y'ix) 

^^/a_ -3/2 _3/2 dvi{s) dv-i{s 



Since vi satisfies the constraint vi < pi we have for x > 0, 

° duijs) ^ dpijs) ^Vaf° ds ^ /g 22) 

.00 {X - Sy {X - S)2 TTt ^\(x - S)2 2t 

where we used the explicit expression ()2.6p for pi. The above inequahty is strict 
since vi cannot be equal to pi on the whole negative real axis. Similarly, we find 

<^^x-3/2, x>0. (3.23) 



,(x-s)2 2(1 -t)" 

Combining (j3.2ip - (j3.23p . we find that (j3.2ip is positive, and therefore (|3.19p is 
strictly convex. It then follows from [46, Theorem IV. 1.10] that the support of 1/2 
consists of one interval, say, [p, q] with < p < q. 

(ii) Suppose that s > is such that i^i = pi and u-^ = p^ on (— s,0). Thus, if x > 0, 



-U^'ix) 



\og\x - y\dpi{y) + j \og\x - y\dpi{y) (3.24) 

-s J —00 



where the second term is real analytic on M^, and for the first term we use the 
explicit expression (|2.5p for pi . Then we obtain for the derivative of (|3.24p (where 
"real analytic" means a real analytic function on R^) 

(-[7^i)'(x) = ^ / + "real analytic" 

+ "real analytic" 





r 


TTt J 




_ ^/a 


r 


TTt J 


—00 


= X 

t 


-1/2 



+ "real analytic" . (3.25) 
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Then by integrating (j3.25p we find that If^^ (x) + is real analytic for x G M"*" 
(including at x = 0). 

Similarly, we find that U"'' (x) + is real analytic on as well. Then V is real 
analytic because of ()2.4p and p.l9p . 

□ 

3.3 Proofs of Theorem 12.11 and Proposition 12.21 

Proof of uniqueness. We rewrite the energy functional £{1/1, 1^2, ^3) as 

Eii^i,U2,i^3) = \l{v2) + \l{y2 - 2vi) + ^I{v2 - 21/3) + j V{x)dv2{x). (3.26) 

It is well-known that if two distinct measures v and /i have finite logarithmic energy and 
J dv = J dfi, then 

I{u - n)>0. (3.27) 

Hence, it follows that the energy functional is strictly convex. The minimizer is 
therefore uniquely determined by the Euler-Lagrange variational conditions listed in 
Proposition [221 □ 

Proof of existence. It is easily seen from (|3.26p ~ ()3.27p that 

E{i^ui^2,i^3) > \l{v2) + j V{x)dy2{x). (3.28) 

Thus the energy functional is bounded from below. As a consequence, we can find a 
sequence {i^i^ni'^2,n,i^3,n) such that 

E{ui n-, 1^2,71, l'3,n) < - + iuf E{ui, 1/2, 1^3), (3.29) 

n 

where both the sequence and the infimum are taken subject to conditions {El)~{E3) of 
the equilibrium problem. In addition, we may assume that 

supp(z^2,n) C\p,q], 

and, in view of (l3TB . (IHTT]) . 

dl'i,n{x) ^ k j = 1 3 

dx ~ ' 

with {} <p < q and K > independent of n. Therefore, there exists a weakly convergent 
subsequence of {yi,n-,T^2,n-,i-'3,n)-, with a certain weak limit {ui, 1^2, 1^3)- By the fact that 
the energy functional E is lower semi-continuous [28J, the vector of measures (z^i, 1^2, 2^3) 
is then the required minimizer of the energy minimization problem. □ 
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Figure 3: The Riemann surface TZ. 



Proofs of Theorem [27T] and Proposition [2721 We already showed that there exists a 
unique minimizer (^ui, /i2, ^3) for the equilibrium problem. The parts (a)-(c) of Theorem 
12.11 and Proposition 12.21 are now immediate by Lemmas I3.1H3.3I □ 

4 Riemann surface and auxiliary functions 

In this section, we shall construct a Riemann surface TZ with the minimizer (/^i, /X2, /U3) 
of the vector equilibrium problem and introduce some auxiliary functions for further use. 

4.1 A four-sheeted Riemann surface TZ 

The Riemann surface TZ has four sheets TZj, j = 1, . . . , 4, which are defined as 

7^l = C \ supp(pi - /ii), 7^2 = C \ (supp(/9i - U supp(^2)), ^^^^ 

7^3 = C \ (supp(/i2) U SUpp(p3 - /is)), 7^4 = C \ SUpp(p3 - /^s). 

These four sheets are connected as follows: TZi is connected to TZ2 via supp(/5i — fii), TZ2 
is connected to TZs via supp{fi2), and TZ3 is connected to TZ4, via supp(/93 — /U3). Each 
connection is in the usual crosswise manner; see Figure O 

The Riemann surface is compactified by adding two points at infinity: one is common 
to the sheets TZi and TZ2, and the other one is common to the sheets TZ3 and TZ4. Since 
the Riemann surface TZ has 4 sheets and 6 simple branch points, it then follows from 
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Hurwitz's theorem (see e.g. [44j) that it has genus 0. This fact will be helpful in our 
future construction of the global parametrix. 



4.2 The ^-functions 

We use Fj to denote the Cauchy transform of the measure ^j, i.e., 

Fj{z) = [ -^dfijix), zgC\ supp(^,), j = 1, 2, 3. (4.2) 

J Z X 

With the above notation, the Sokhotski-Plemelj formulas now read 

F,- +(z) + Fj,.{z) = 2PV / -^dfxjix), (4.3) 

J Z X 

Fj,^{z)-F,,4z) = -2m^, (4.4) 

for z G supp(^j), where PV denotes the Cauchy principal value. The ^-functions are 
then defined as follows. 

Theorem 4.1. The function : Uj=i ~^ defined by 



-^i(^) + ^ + 7(T^' ^e7^l, 
Fi(z)-F2(z)-^ + ^, zG7^2, 
F2iz)-Fs{z) + j^^, ^G7^3, 



(4.5) 



has an extension to a meromorphic function (also denoted by S,) on TZ. In ()4.5p . the 
principal branch of the square root function ^/z is used, i.e., with a branch cut along 
(-00,0]. 

Proof. We first need to show the analyticity of on each individual sheet TZj. To see 
this, we note that = pi on M~ \ supp(pi — /^i) = (— ri, 0). It then follows from (j4.4p 
and the definition ()2.6p of pi that 



dx dx 

for X G (— '"I, 0). By the definition of ^, this implies that ^ has an analytic continuation 
to (— ri, 0) on both IZi and 7^2- Similarly, we have by ()4.4p and ()2.8p 

F„(.)-F,_(.) = ^(-^--^). (47) 
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for X G (— r2,0), which gives the analyticity of ^ on the sheets TZ^ and TZ^. 

It remains to check the analyticity of ^ if one moves from one sheet to another via 
a specific cut. This is always the case in view of the variational conditions for the 
equilibrium problem. For convenience, we will only show that is analytic when one 
crosses supp(/?i — ^i), since the other cases can be verified similarly. 

We start with differentiating U^^ on supp(/9i — fii) and obtain 

2—U^'^{x) = -2PV / dfxiis) = - (4.8) 

ax J X — s 



where the second equality follows from Sokhotski-Plemelj formula (|4.3p . On the other 
hand, differentiating the right-hand side of the variational condition (|2.17p . we have 

2i-[/m(^) = ^c//^2(^) = _ f = -F2{x). (4.9) 

ax ax J X — s 

Combining (|4.8|) and (|4.9|) leads to 

Fi,+ (x) + Fi,_(x) = F2(x), X e supp(pi - A^i). (4.10) 

This, together with the fact that y/x_ = — -v/x+ for x G M~, gives the analyticity of ^ 
when one crosses supp(pi — /ii) and moves from IZi to 7^2- D 

As a consequence of the above theorem, the following result about the asymptotic 
behavior of Fj[z) at infinity is easy to prove. 

Lemma 4.2. There exist two constants ci and C3 such that 

Fii^) = iz-^ + 0(^]. (4.11) 



1 


Cl 


Yz 


^3/2 


1 




- 4 




z 




1 


C3 


2I 


z3/2 



F2(z) = - + 0(^), (4.12) 
^3(.) = --T^ + 0f^\ (4.13) 



as z 00. 



We use ^j{z) to denote the restriction of ^ to TZj. The relations between and the 
solution to the equilibrium problem can be established in the following way. 

Proposition 4.3. We have that 

dfJLi{x) = ^XiiAx)-iiAx))dx + dpi{x), xGR-, (4.14) 

ZTll 

dfi2{x) = ^iC2,+ {x) - 6,-(x))t^x, X E M+, (4.15) 

dfi3{x) = A^{^3^+{x)-C3,^{x))dx + dp3{x), xeR-. (4.16) 

2m 
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Proof. This follows by a combination of (|4.5|) . the Sokhotski-Plemelj formula (|4.4|) . and 
the definitions ()2.6p and (|2.8p of pi and p2- □ 



The Riemann surface has two points at infinity and the meromorphic function (j4.5p 
takes on two distinct finite values at these points (namely tp^rjj and 0). Thus ^ is not 
a constant, and being a meromorphic function on a compact Riemann surface it must 
have at least one pole. Prom ()4.5p it is clear that the only possible poles are at the values 
of z = that correspond to branch points. It thus follows that at least one of ri,p, and 
r3 is zero, and we can identify the various cases in terms of the structure of TZ. 

In Case I we have ri = = and p > 0. It is only in this case that ^ has two poles. 
In all other cases, ^ has exactly one pole, and then ^ is actually a conformal map from 
TZ onto the Riemann sphere. 

4.3 The spectral curve 

From Theorem 14.11 it follows that the four function (,j{z), j = 1,2,3,4 are solutions 
to an algebraic equation of degree four 

= p(e, z) = {^- ei)(^ - 6)(e - - u) 

= + A{z)^^ + B{z)f + C{z)C + D{z), (4.17) 

with coefficients A(z), B(z), C{z), D{z) that can only have a pole in z = 0. In fact, 
because of the behavior of (|4.5p at z = 0, we have that A(z) has no pole at z = (and 
therefore is a constant), B{z) and C{z) have at most a simple pole at z = 0. Also D{z) 
has at most a simple pole at z = 0, unless we are in Case I, in which case D{z) has a 
double pole at z = 0. Indeed, in this last case we have that all functions behave like 
as z — )• 0, and then D{z) = ^i{z)^,2{z)£,3{z)£,4{z) behaves like z~'^ as z — )• 0. 
Straightforward calculations using ()4.5p and Lemma 14.21 vield in all cases 

A{z) = -6 - 6 - 6 - ^4 = -j^rzT)^ (^-^^^ 
B{z) = 66 + 66 + 66 + 66 + 66 + 66 

1 b a 1 , , ^ 



t2(l-t)2 z{l-tf zt^ zt{l-t) 

^(^) = -;^^(rri)4 + iW^' ^'-''^ 

where c = in all cases but in Case I. 

Lemma 4.4. The spectral curve is given by ()4.17p with coefficients given by (j4.18p - 
()4.2ip . The constant c that appears in (|4.2ip is equal to 

-i)' ^nCaseI, ^^^2) 
in all other cases. 
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In Case I, we have 

^=(l-t)^ + tVb-^/2t{l-t)>0, ^^^^^ 
^ = (1 - t)^/a + t\/6 + ^2^(1^^. 

Proof. We already saw that the spectral curve is given by ()4.17|) . ()4.18p - ()4.21|l and that 
c = in all cases, except in Case I. 

Now suppose we are in Case I, so that ri = = and p > 0. Consider the following 
functions 

[z£,i{z''), Re2;<0, [2^2(2), Rez<0, 

jz^siz^), Rez>0, izU{z^), Rez>0, 

[2,^4(^2), Rez<0, [z£,3{z'^), Re2;<0. 

Since, we are in Case I, we have that all functions rjj are analytic across the imaginary 
axis. Also the singularities at z = are removable because of the extra factors z in 
()4.24p and ()4.25p . It follows that rji and % are analytic in C \ [^/p, ^/q], while r]2 and 774 
are analytic in C \ [— y^, —^/p]- 

In fact the two functions r/i and 773 give a meromorphic function on the two-sheeted 
Riemann surface with cut along [y/p, ^/q], whose only pole is at infinity on the first sheet. 
Since 

as z — )• 00, it then follows that 771(2;) and 7/3(2) are the two solutions of the quadratic 
equation 

Similarly, 7/2(2) and 7/4(2) are the two solutions of 
r? - ( ' + ^ - 7/ + ( ^ ^ + -^—] = 

' l^t(i-t) t i-tj'^y tii-ty t{i-t)^ 2t{i-t)j 

In particular, it follows that 



/ ^ / N / N / N ^z^ {Vab - i 

rii{z)m{z)m{z)V4{z) = --jTTz -77 + 



1\2 



t2(l-t)4 t2(l_i)2 
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Then in view of (|4.24p and (|4.25|) we obtain 



2 



t2(l -t)4 t2(l _i)2 • 

Comparing this with ()4.2ip we indeed find that c is given by ()4.22p in Case I. 

The branch points of (j4.26p are equal to ^/p and y^. Simple calculations then lead 
to (1123]). □ 

Remark 1. We notice that the spectral curve in Case I is equivalent with the one in 
jlHl , while in Cases II and III it is equivalent with the one in flS^J ( with ^/a, y/h instead 
ofa,b). This implies that the limiting distribution of the non-intersecting squared Bessel 
paths is obtained from the one of the corresponding non-intersecting Brownian motions 
by a simple squaring transform. 

One could extend the above ideas even further by also lifting the vector equilibrium 
problem into its square-root form. The lifted measures are then supported on the real 
or imaginary axis. This approach can be used to obtain an equilibrium problem inter- 
pretation for the spectral curve in 1 13^ . which allows in turn to obtain an alternative 
derivation of the results in that paper. We will not go into the details. 

Proposition 4.5. We are in Case I if and only if ab > 1/4. 

Proof. The spectral curve depends continuously on the parameters a, b and t. Therefore 
the change from Case I to one of the other cases can only occur when ab = 1 /4. 

If we are in Case I, then ^/p is given by (j4.23p which in particular implies that 
(1 - t)y/a + tVb - y/2t{l - t) > 0. For t = 1/2 and o = & < 1/2 this inequality is not 
satisfied and so for these values of parameters we are not in Case I. Then it follows that 
for all a and b with ab < 1/4 we are not in Case I. 

If we are not in Case I, then c = and the spectral curve (|4.17p has four real branch 
points — ri, — r2, p and q. For the specific values t = 1/2, a = 1 and b = 1/2, we 
calculated with Maple that there are two non-real branch points. Thus for those values 
of parameters we cannot have c = 0, and so are in Case I. By continuity we are in Case 
I for all a and b with ab > 1/4. This completes the proof of the proposition. □ 



4.4 Phase transitions 

Let us assume that the ending point b is fixed. In figure H) we give a phase diagram 
divided into three different regions according with the limiting distribution of the non- 
intersecting squared Bessel paths as n — )• oo. The horizontal axis of the figure denotes 
the time t G (0, 1) and the vertical axis denotes the starting point a. 

In the case ab < 1/4, for t small, we are in Case II. The paths start at the positive 
value a at t = and they initially remain positive. The limiting mean density is sup- 
ported in the interval [p,q] with p > (see Theorems 12.11 and 12. 3p . For this case there 
is only one active constraint. 
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Figure 4: The left figure shows the phase diagram for the non-intersecting squared Bessel 
paths for b = 1/4. The sohd curve is given in ()4.27p . The right figure shows a simulation 
of 50 non- intersecting paths for a = 1/3 which corresponds to the dashed line in the left 
figure. 



We have computed the discriminant Disc{z) of ()4.17p with coefficients given by 
()4.18p - ()4.2ip and c = in ()4.2ip with respect to ^ using Maple and we obtain 

Disc{z) = a4z'^ + a3Z^ + 02^^ + 4(4a6 - l)(a - 2at + at^ + ht^ + - tfz 

for some constants 04, 03 and 02- We do not include the explicit expressions of 04, 03 
and a2 because they are too long and not too interesting. Therefore there are exactly 
four zeros Disc{z) which correspond to the branch points of the Riemann surface. 

For each ab < 1/4 there are two critical times ti^2 where Disc{z) has a double zero 
at z = 0, namely 

2a + l±Vr^4^ 
'''' - 2(a + 6+l) • 

The critical times ti and ^2 are the only two values of t for which the coefficient of z in 
Disc{z) is equal to zero. The time ti corresponds to the time where the first paths touch 
the wall. If ti < t < t2 then we are in Case III. Here the limiting mean density of the 
paths is supported on the interval [0,q] (see Theorems 12.11 and 12. 3p . Both constraints 
are active in this case. If t > t2 then we are back to the Case II. 

The figure on the right of Figure U] shows the transitions between Case II and Case 
III for the choice of the parameters a = 1/3 and 5 = 1/4. This choice corresponds to the 
dashed line in the left figure. The critical times ti and ^2 where the transition occurs 
are given by (|4.27p . 

On the other hand, if a5 > 1/4 then the paths start at the positive value a at t = 
and they remain positive. In this case the limiting mean density is supported in the 
interval [p, q] with p > 0. No constraint is active in this case. 
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4.5 The A-functions 

For the measures j = 1,2,3 that minimize the energy functional (j2.3p . we first 
define the associated g-functions as 



gj{z) = J \og{z-s)dtij{s), (4.28) 

where we take the principal branch of the logarithm, i.e., 

log(2; — s) = log 1^ — s| + i arg(2; — s), arg(2; — s) G (— vr, vr). (4.29) 

On account of Theorem 12. H it is easily seen that gi and 53 are analytic in C \ , while 
g2 is analytic in C \ (— oo,g]. 
Noting that 

g'^{z) = F,{z), (4.30) 
we find the asymptotic behavior of gj (z) at infinity as given in the following lemma. 
Lemma 4.6. With the constants ci and C3 given in Lemma \4-S\ we have that 



9i{z) 


= ^log 


2ci 


92{z) 


= log z 


Hi 


93{z) 


= - log 
2 ^ 


^ 2C3 



(4.31) 
(4.32) 



+ O ( i ) , (4.33) 



as z —)• 00. 



Proof. The asymptotics of (72 follows from the fact that ^2 is a probability measure with 
compact support. 

The formulas ()4.3ip and ()4.33p for gi and g^ can be derived by integrating the 
asymptotic formulas for Fi and given in Lemma 14.21 The integration constant is 
zero, due to the fact that 

j log{z - s)dfij{s) = log(z)||^j|| + j log (^1 - dfij{s), j = 1,3, 

and / log(l — s/z)dfij{s) — )• as z — )• 00. □ 
The A-functions are then defined as the following anti-derivatives of the ^-functions: 

'^^ + 777^-'. (4.34) 



Ai(z) 


= -gi{- 


Hz) 


= gi{z) 


Hz) 


= 92{z) 




= 93iz) 



t t{l-t) 



t t{l-t) 
2\/te 



2Vbz 



(4.36) 
(4.37) 
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where / G M is the variational constant in Proposition 12.21 Note that Ai and A4 are 
analytic in C \ (— oo,0] while A2 and A3 are analytic in C \ (—00,(7]. Some properties of 
A-functions are shown in the following two lemmas. 

Lemma 4.7. We have 

Ai,+ (x) = Ai xG(-ri,0), 

A2,+ (x) = A2 -(x) - 2TTi, X £ (-ri,p), 

h,+ {x) = X3-{x) + 27ri, xG(-r3,p), 

A4,+ (x) = A4-(x), xG(-r3,0), 

^i,±ix) = A2,=f(x) =f vri, X £ Ai, 

M,±ix) = A3,=p(2;), X £ A2, 

A3,±(a;) = A4,q;(x) =F TTi, X £ A3, 

and 

Re(Ai(x) - A2(x)) < 0, xG(-ri,0), 

Re(A2(x) - A3(x)) > 0, x G M+ \ A2, 

Re(A3(x)-A4(x)) <0, xG(-r3,0). 

Proof. These are reformulations of the variational conditions of the equilibrium problem 
in Proposition 12.21 □ 

Lemma 4.8. The X-functions have the following asymptotics for z ^ 00: 



X / \ ^ 2Jaz , 1 ^ 2ci ^, 

^^(^) = ^(r37) + ^-^-2^°^^-7i + ^(^ ' 

, / ^ 2v/ai , 1 , 2ci ^, 

^^(^) = ^(i^-^-^-2^°^^ + 7i + ''^' ' 

2^z \ 2c3 ^ _i, 

A3(.) = y3^ + 2log^-;7^ + 0(- ), 



Proof. Obvious from the definitions (fi3i]) - (f071) and from (fOT]) - p:33|) . □ 

4.6 The </)-functions 

We end this section with the introduction of the i?i>-functions, which will be used to 
simplify the jump matrices during our steepest descent analysis. Recall that is the 
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restriction of ^ to the sheet TZj of the Riemann surface TZ. We define 

<t>i{z) = \ j\ii{y)-Uy))dy. -2GC\((-cx),-ri]U[p,oo)), (4.38) 

02(^) = \j^ iUy) - Uy))dy, ^ e C \ (-00, q], (4.39) 
Mz) = l r {Uy)-Uy))dy, zeC\{{-oo,-rs]u[p,oo)). (4.40) 

5 Steepest descent analysis for Y{z) 

In this section, we will perform the steepest descent analysis of the RH problem for 
Y. It consists of a series of explicit and invertible transformations. In Section El these 
transformations will be used to prove Theorem 12.31 

To this end, we mainly follow the theme laid out in ^37j, although the argument is 
somehow more involved since our RH problem is of size 4x4 whereas the dimension 
treated in [37] is 3 x 3. Furthermore, we have more generic cases to consider here. 



5.1 First transformation Y X 

It is the aim of this transformation to simplify the block matrix ()2.30p appearing in 
the jump condition ()2.29p for Y. The idea is to use the special properties of modified 
Bessel functions, in a very similar way as in [37j. To start with, we set four functions: 

y,{z) = z("+l)/2/„+i(2Vi), y2{z) = z(°+l)/2K„+i(2Vi), 

ysiz) = z-("-i)/2/,_i(2Vi), y^z) = z-("-i)/2^„_i(2Vi), 

where Ki, is the modified Bessel function of the second kind; see [U Section 9.6] for its 
main properties. It is readily seen that y3{z) is well defined in the whole complex plane, 
while the function yj{z),j = 1,2,4 is analytic in the complex plane with a branch cut 
along the negative real axis. Indeed, it follows from formulas 9.6.30-9.6.31 in [1] that 
their jumps on M" are given in the following ways 

2/i,+ (x)=e2--yi,_(x), 

2/2,+ {x) = y2A^) + ivre^"-yi,_ (x) , (5.2) 
2/4,+ (x) = e-2^"^2/4,-(a;) +ivre-^°^y3(x), 

for X < 0. In view of the derivative properties of the modified Bessel functions JV, 9.6.26], 
we also have 

y[{z) = z-/^U2^z), y',{z) = -z"'^ K^{2^z), 
y'.iz) = z""/2/„(2Vi), y'Az) = -z-'^/^K^il^z). 
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A combination of (|5.1|) . (j5.3|) and (|2.30|) implies that we can represent the weight 
matrix W(x) in terms of the functions yi, 1/3 and their derivatives as 



2n\/ax 



= '""^ r A) '■f-'!/3(r|.)) . (5.4) 

where 

n = ri_„ = — T3 = r3,„ = J— ^. (5.5) 
For later use, we also need the following Wronskian relations: 

y[{z)y2{z) - yi{z)y'^{z) = /2, y'3{z)yi{z) - y3{z)y'^{z) = z-" /2, (5.6) 

for z G C \ MT] see formula 9.6.15 of |T]. 

The first transformation Y ^ X is then defined by (compare with j37j ) 



where 



X{z) = CiY{z)dmg{Ai{z),A2{z)), (5.7) 



A,{z) = rr"z-"/2 ^ ; \^ , (5.8) 

-i,rriy2(r^) T^^yiirlz)^ 

-ny'^ir^z) TTiTsy'^irlz) 



and with Ci the constant 4x4 matrix 



Ci = diag V2vrri 



1 



4("+l)"-l jl'V2^1^ , 4(a- 1)2-1 
16ri / \ 16t3 



Note that, on account of (15.61). we have 



detAi = —i, detA2 = 2TTTsi. (5.10) 

2ttti 

Hence, it is easily seen that detX = 1. Now, we have 

Proposition 5.1. The matrix valued function X{z) defined by ()5.7p is the unique solu- 
tion of the following RH problem. 
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(1) X{z) is analytic in C \ M. 

(2) For X G M, ^i^) satisfies the following jump condition 



X+{x) = X_{x) 



I + e *(!-*)£': 



diag(e" 



23; 



21 



43) 



ifx>0, 
ifx<0, 



(5.11) 



where we denote by Eij the 4x4 elementary matrix whose entries are all 0, except 
for the {i,j)-th entry, which is 1. 



(3) z oo, z G C \ M, we have 
X{z) = {l + o{^^^dmg 



1 




X diag ( ^"/2g-2nv^A ^ ^n/2g2nv^/t^ ^-n/2g-2nV^/(l-t) ^ ^-n/2g2nv^/(l-i) ' 



XZ^{x)X+{x) 



Yl^{x)Y+{x) 



uniformly valid for z hounded away from the negative real axis. 
(4) X{z) has the following behavior near the origin: 

X(z)diag(|z|"/Mz|-°/Mz|"/Mz|-°/2)) = 0(1), if a > 0, 

X(z)diag((log|z|)-i,l,(log|z|)-i,l) = 0(1), i/a = 0, 

X(z) = 0(z"/2), X-i(z) = 0(z"/2)^ i/-l<a<0. 

Proof. The proof follows [37j. It is easily seen from ()5.7p ~ ()5.9p and ()2.29p that 

^A^lix) 
A^^Jx) 

A'{}_{x)W{x)A2,+ {x) 
h 
A^^_{x)Ai^+{x) 

A^^{x)A2,+ {x) 

In view of (j5.4p and (|5.6p . we obtain from a direct calculation that 
A-{'l{x)W{x)A2,+ {x) 

r^'^v'xijlx) 




if X > 0, 



if X < 0. 



e (i-*)*^ii(x) 



e (!-*)* 




vAtIx) 



-a/2 




(5.12) 



(5.13) 



(5.14) 
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for X > 0. Also, by (|5.2p and i\5.'6\i . it can be shown that 



-1 



and ^2 i(x)^2,+ (a;) 



-1 



, (5.15) 



for X < 0. The jump condition ()5.1ip now follows from a combination of ()5.13p - ()5.15p . 

The large z behavior of X(z) shown in item (3), follows very similarly as in [37J. The 
relevant results are 



A,iz) 



1 







' ' 16ri ^ 



+ 0{z-') 



-^^3/4 / 1 



V2 



i 1 



(5.16) 



and 



A2iz) = 



' 4(a-l)^-l 
16r3 







+ Oiz-') 



x-<^3/4 / 1 



V2 



i 1 



(5.17) 



as 2 — )• oo, uniformly for z bounded away from the negative real axis. Substituting the 
asymptotic formulas (j5.16p and (j5.17p of Ai(z) and ^2(2) into (j5.7p . the asymptotic 
behavior of X{z) at infinity is immediate. 

Finally, the known behavior of the modified Bessel functions near zero given by 
formulas 9.6.7-9.6.9 in [T] yields 



1 



1 



r(a + 2)' 



r(a + 1) 
-ir(a). 



y2{z) ~ -r(a + 1) 



Q > 0, 

a = 0, 



^log(z), 

-ir(-Q)z°, Q < 0, 



and 



^3(2) 



T{ay 



ir(a-l)zi-", a>l, 

2/4(2) ~ { _1 log(2;), a = 1, 

|r(l — a), a < 1, 



2/3 



2/4(2) 



1 



r(a + i)^ 



-ir(a)z-", a > 0, 
ilog(z), a = 0, 
-ir(-a), a<0. 



The behavior at zero of X(z) in item (4) then follows from a straightforward calculation. 
This completes the proof of Proposition 15.11 □ 
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5.2 Second transformation X ^ U 

The second transformation X ^ U \s defined by 



U{z) = C2X{z) diag \^'V''> W=T)J ^ ^^-^yi-W^T)) ^ gnA3(^j^ gnA4(^) J ^ (5 jg) 

where C2 is the constant matrix 

C2 = (/ + 2inci^2i + 2inc3Eu) diag(e"', e"', 1, 1), (5.19) 



with ci and C3 as in Lemma 14.21 and with / the variational constant in Proposition 12.21 
In view of the jump relations (jS.lip for X and the definition ()5.18p . we see that the 
jump relations for U are given by 

f A-'^ (I + E23)A+, onM+, 
U+ = U-{ - ^ ' (5.20) 

[ AI^ (diag(e-*™,e*™,e-*™,e*™) -^21 -£^43) A+, onM", 

where we define the diagonal matrix 

A(z) = diag(e"^i(^\e'"^2(^)^gnA3(^)^gnA4(^))_ (5^21) 

The jump matrix for U can be reformulated with the help of the ^-functions in 
()4.38p - (|4.40p as stated in the following proposition. 

Proposition 5.2. The matrix valued function U (z) defined by (j5.18p is the unique 
solution of the following RH problem. 

(1) U{z) is analytic in C \ R. 

(2) For X gW, we have that [/+ = U^Ju, where 

on A 

1 e2"*2>- / ' / 
Ju = I + e~2"02,+ on M+ \ (0 U A^) , 

Ju = diag \ \ 



Ju=\\ ° if'' ° I I on Ai \ A 

g«7rag2n(/)i,„ / I _g2n03 gjvra 

Ju = diag(e-*™, e^™, e"^™, e*™) - e^'^'*'^E2i - e^'^^^E^^ on \ A 






with Aj, j = 1,2,3 being defined in (j2.12p and Aj denotes the closure of A 
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J- 



(3) As z ^ oo, z G C \ M, we have that 



U{z) = (^/ + oQ^^ diag(z-l/^zV4^^-l/4^^l/4 



IW .. / MA ma\ I /I M /I ^ 

^ " VV V ' V 1. 

(5.27) 

uniformly valid for z hounded away from the negative real axis. 

(4) U{z) has the same behavior as X{z) near the origin; see ()5.12p . 

Before proceeding to the proof of Proposition 15.21 we remind the reader that, due to 
our assumption ()2.13p . the jump condition ()5.25p only appears in Cases II and III, while 
(j5.26p only appears in Case III. 

Proof. The jump conditions in item (2) are obtained in a straightforward way from 
Lemma 14.71 and the definitions of i;^- functions in ()4.38P " ()4.40p . For convenience, we give 
a proof of (j5.22p . the other jumps can be proved in similar ways. 
If G A2 , it follows from ()5.20p and Lemma 14.71 that 

/ /gn(A2,+-A2,-)(x) I \ \ 

Jt/(x) = diag 1, . 'M- (^-2^) 

On account of our definition of A2,3 in (j4.35p and (j4.36p . it suffices to show that 

g2,+ {x) - 52,-(x) = -202,+ (x) = 2(t>2,-{x), (5.29) 
for X G A2. To see this, we obtain from (|4.28p and ()4.15p that 

ff2,+ (a:) -52,-(x) = 27ri rd/X2(s) = / (6,-(s) - 6,+ (s))ds. (5.30) 

J X J q 

The relation (|5.29p now is immediate, with the aid of ()5.30p and the fact that ^2,- {s) = 
provided s G A2. 

To prove the asymptotic behavior of U near infinity, we note that the asymptotics 
of X and the large z behavior of the A-functions given in Lemma 14.81 yield 

X{z) diag I e 



(5.31) 

as z — )• 00, where 

. dia, (1 - ^) , (, , H) , 1 - 1 , 51) , 0(i) . 
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where the matrix denoted with O(^) is diagonal as well. 

Now we can move C{z) in ()5.3ip to the left as in the proof of Proposition 15.11 The 
result is that (|5.3ip is equal to 

Then ()5.27p is immediate by the definition of U in ()5.18p . 

The behavior of U near the origin follows from the corresponding behavior of X, and 
the fact that the A-functions are bounded at the origin. □ 



5.3 Third transformation U ^ T 

Our next transformation involves factorizations of jump matrices and the so-called 
lens opening, which is a key ingredient of the steepest descent analysis. The main goal 
of this step is to convert the highly oscillatory jump matrices into a more convenient 
form (e.g., constant matrices) on the original contour, thereby introducing extra jump 
matrices with exponentially small off-diagonal entries on the new contours. 

The existence of lenses is guaranteed by the next lemma. 

Lemma 5.3. (Lenses around Aj.j There exist Jordan curves A^, Ag connecting p to q 
in the upper and lower half plane, respectively, such that Re (^2 (-2) < for every non-real 
z in the region enclosed by A^, A2 . This region is called the lens around A2. 

Similarly, for j = 1,3 there exist Jordan curves A^, Aj connecting —00 to —rj in 
the upper and lower half plane, respectively, such that Re (jyjiz) > for every non-real z 

in the region enclosed by A^, Aj , j = 1,3. This region is called the lens around Aj. If 

± 
1 • 



we are in Case I or r^ = ri in Case III, then we may (and do) assume that A^ = A^ 



Proof. We first prove the claim for A2. Recalling that S^2,+ ix) = S,3-{x) for x £ A2, it 
follows from ()4.15p and the definition ()4.39p of (j)2{x) that 

/•X 

(j)2,+{x) = -(j)2,^{x) = Tri / dfi2, X e A2. (5.32) 

Jq 

Hence, 4)2,± is purely imaginary on A2. Furthermore, by differentiating the above formula 
with respect to x on both sides, we obtain 

-^Im 02 +(2;) > and -^Im 02 -(x) < 0. (5.33) 
dx ' dx ' 

Therefore, an appeal to the Cauchy-Riemann equations yields the existence of a region 
around A2, such that Re 02 (-2) < for every z ^ {p, q} in that region. 

To prove the claim for Ai, note that ^i^+(x) = ^2,- (3;) for x G Ai. Then we obtain 
from (OHj) and (liTijl that 

/X 
d{fii-pi), X e Ai. (5.34) 
-ri 
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In view of the constraint condition fii < pi, it is immediate that 

-^Im (/)i,+(x) < and 4-1™ 0i > 0. (5.35) 

Again, by the Cauchy-Riemann equations, our claim holds. The claim for A3 can be 
proved similarly. □ 

An illustration of Lemma 15.31 is shown in Figure [5l 

The 2x2 middle block of the jump matrix (|5.22p on A2 admits the factorization 



e2n</.2,+ 1 \ / 1 




(5.36) 



Similarly, for the 2x2 upper left and lower right blocks of the jump matrix ()5.24p on 
Ai and A3 we have the factorizations 

j = l,3. (5.37) 

Based on the factorizations (|5.36p - ()5.37p . we define 

T{z) = U{z) (l ^ e2"'^2(^)E32) , (5.38) 
for z in the domain bounded by A^ and A2; 

T{z) = U{z) (l ± e^^^e-^^-^i^^^^^ia) , (5.39) 
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for z in the domain bounded by A]*^, A3 and Ai; 

T{z) = U{z) (^I ± e^^™e-2"<^i(^)^i2 ± e±^™e-2"*3(^)_E34^ , (5.40) 

for z in the domain bounded by A^ and A3, and we let 

T{z) = U{z), (5.41) 

for z outside of the lenses. 

It then follows from straightforward calculations that T(z) is characterized by the 
following RH problem: 

Proposition 5.4. The matrix valued function T{z) is the unique solution of the following 
RH problem. 

(1) T is analytic m C \ (E U A]^ U Af U A^). 

(2) For each of the oriented contours shown in Figure\^ T has a jump T+ = T^Jt, 
where 

( (0 A 

Jt = diag 1, , 1 , on A2, 

V v-1 V J _ 

Jrp = 1 + 6-^''^^'+ E23, on M+ \ (0 U A2), 

Jt = diag(e-^™,e*™,e-*™,e^™) - e'^^''^'E2i - e2"<^3^43, on (-ri,0) {Case III), 

,0 l\ / e-^'^" \ \ 

Jt = diag , , on Ai \ A3 [Cases II and III), 

-1 Oy y-e2"<^3 ^inaj J 

Jt = diag | | , , on A3, 

-1 0) \-l o)J 

Jr = / + e2"'^2^32, onAf, 

Jt = I- e^*"^e-2"<^i^i2, on Af , 

Jt = I- e^*"^e-2"'^3E34, on Af , 

Jt = I- e^*"^e-2"<^iEi2 - e^^°^e-2"<^3^34, on Af = Af {Case I or r^ = ri in Case 

(3) As z ^ 00, we have 

T{z) = (^/ + oQ^^ diag(z-l/^zl/^z-V^zV4^ i=diag 
uniformly for z G C \ M. 

34 



i 1/ ' U 1 



(4) T{z) has the same behavior near the origin as U{z) (and X{z)), see (|5.12|) . provided 
that z — )• outside the lenses that end in 0. [The behavior inside the lenses is 
different but will not be needed for us.] 

(5) T is bounded at the finite endpoints of Aj, j = 1,2,3 other than the origin. 

Proof. All the properties follow from straightforward calculations. 

To show that the asymptotic behavior of T given in item (3) holds uniformly up to 
the negative real axis, one needs to trace back the transformations Y X U i— t- T. 
It turns out the entries of T are actually linear combinations of the modified Bessel 
functions, whose asymptotic expansions are uniformly valid up to the negative real axis. 
We omit the details here. □ 

5.4 Global parametrix 

By the above constructions, it is easily seen from Lemmas 14.81 and 15.31 that the jump 
matrices Jt in the RH problem for T{z) all tend to the identity matrix exponentially fast 
as n — 7- oo, except for the jump matrices on Aj, j = 1, 2, 3. Hence, we may approximate 
T by a solution of the following 4x4 model RH problem: 

(1) Na is analytic in C \ (A^ U A^ U A^). 

(2) Na has continuous boundary values on Aj, j = 1, 2, 3 and satisfies A^a,+ = Na^-J]\j^ , 
where 

on A3, 

on Ai \ A^ (Cases II and HI), 

on (-ri,0) (Case HI), 
on A2. 



(3) As z ^ 00, z G C\M-, 

No.{z) = + diag(.-V4,,i/4^^-i/4^^i/4^ _l^^i^g 

We construct the solution to this RH problem with the help of meromorphic dif- 
ferentials, following the idea in |40) . We will give the construction for Case HI; the 
modifications for Cases I and II will be briefly commented at the end of this section. 





35 



Consider the Riemann surface in Figured Denote with p, q (sheets 2 and 3), — ri 
(sheets 1 and 2), —r^, (sheets 3 and 4) the finite branch points. Also denote with 001^2 
the point at infinity that is common to the first and second sheet, with 003^4 the point at 
infinity that is common to the third and fourth sheet, and with Oi and O4 the origin on 
the first and fourth sheet respectively. Recall that we are working in Case III, so p = 
and min(ri,r3) > 0. 

Let Wj, i = 1,2,3,4 be four meromorphic differentials (of the third kind) on the 
Riemann surface with simple poles at the above listed points (and nowhere else), and 
with residues at these points given by the following table: 



X 


P 


q 


-n 




001,2 


003,4 


Oi 


O4 


Res(a;i, x) 


-1/2 


-1/2 


-1/2 


-1/2 


1/2 


3/2 


-a/2 


a/2 


Res(w2, 2;) 


-1/2 


-1/2 


-1/2 


-1/2 


-1/2 


5/2 


-a/2 


a/2 


Res(a;3, j;) 


-1/2 


-1/2 


-1/2 


-1/2 


3/2 


1/2 


-a/2 


a/2 


Res(a;4, x) 


-1/2 


-1/2 


-1/2 


-1/2 


5/2 


-1/2 


-a/2 


a/2 



The differences among Wj are that their residues at 001,2 and 003,4 are all different. Such 
meromorphic differentials exist because the sum of the residues is zero. The fact that 
our Riemann surface has genus implies that each meromorphic differential is unique; 

cf. mM- 

We then build the first row of with the aid of meromorphic differential coi intro- 
duced in (|5.42p . Let pi be an arbitrary reference point on the first sheet of the Riemann 
surface, chosen on the positive real line such that pi > q. Define four functions 

Uj{z) = wi, i = 1, . . . ,4, 

where we denote with zj the point z lying on the jth sheet. The integration path from 
pi to Zj obeys the following rules: It can only move from one sheet to another by going 
from the lower side of Aj on sheet j to the upper side of Aj on sheet j + 1- Moreover, 
the path is not allowed to go through a pole of loi, and it must not intersect the intervals 
[— ri,0] on the first and second sheets, nor the intervals [-^3,0] on the third and fourth 
sheets, except possibly at the endpoint of the path. 

An illustration of the integration path for U3{z) is shown in Figure [6l 

With the above conventions, the functions Uj{z) satisfy the following jump properties 

ui,-l- = ui,_ + 27ri Res(a;i, Oi) = ui,_ — ani, 
U2,+ = ii2,- — 27ri (Res(a;i, Oi) + Res(a;i, 001,2) + Res(a;i, 
= U2,- + avri, 

U3,+ = U3- — 27ri (Res(a;i,04) + Res(a;i, 003,4) + Res(a;i, 

= U3- - 27ri (1 + a/2) , 
U4,+ = ii4,_ + 27ri Res(a;i, O4) = U4- + ani, 



-n)) 



-rs)) 



on (-ri,0), 

on (-ri,0), 

on (-r3,0), 
on (-r3,0). 



(5.43) 
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Figure 6: Integration path for U3{z). 



U2,+ = 


Ml _, 


U2,~ 


= ui,+ 


+ 2Tri Res(cL!i, 001^2) = ^^1,+ + ^ri, 


on Ai 


U3,+ = 




Us- 


= U2,+ 


— 27ri Res(a;i, q) = n2,+ + vri, 


on A2 


= 


U3-, 


ui - 


= U3,+ 


+ 27ri Res(a;i, 003^4) = n3^+ + Sni, 


on A3 


Now set 


















Vj{z) -- 


= exp{uj{z)), j = l,...,4, 





Each of these jumps can be easily verified by deforming the integration path of uj^^—uj^^ 
to a closed loop and checking which poles are inside (or outside) the loop. Similarly, we 
have 



(5.44) 



(5.45) 

and consider the row vector v{z) := {vi{z),V2{z),vs{z),V4{z)). By ()5.42P " (|5.44p we find 
that 

v+=v_JAr^, onAiUA2UA3, (5.46) 

in Case III. Thus v(z) satisfies the required jumps for the row vectors of Na- 

Next we discuss the asymptotics of v(z) for z — t- 00. Taking into account that c«i^2 
is a branch point on the Riemann surface, its local coordinate on the Riemann surface 
can be chosen as w{z) = l/\/z. Therefore 

Uj{z) = Res(cLii, 001^2) log w(z) +7j + 0{w{z)) 
= -^logz + 7j +0(z"^/^), 2; -^00, 
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for j = 1,2, where 71 and 72 are certain constants. From the relation U2,+ = ui- on 
(—00, — ri) we find that 72 = 71 + Res(ti;i, 001^2)71"^ = 71 + ni/2. Similarly we also have 

Uj{z) = --log z + jj + Oiz~^^'^), z — > 00, 

for j = 3, 4, where 73 and 74 are certain constants. Combining all of this and using 
and I^Mh . we find that 



00. 



The behavior of v{z) at the origin can be found from similar arguments as above and is 
given by 

v(z) = O {z~^/\ z-^l\ z-^l\ z^/') , z^O. 
Finally, we also find that 

^,{z)=0{{\z-K\)-^/^), asz^K, (5.47) 

if K is a branch point other than the origin. 

Summarizing, we see that ^^^^(z) satisfies all the constraints for the first row of 
Nq,. The other meromorphic differentials coi, i = 2,3,4 in Table [532] can be used to 
build the i-th row of in a similar way, we omit the details here. 

Our construction of A^^^^ leads to 

Naiz) = 0{{\z- k\)-^/^), as z^k, (5.48) 

if K is a branch point other than the origin, and 

Na{z) diag(z"/2^ z'^/^, z-"/^) = 0(1)^ as z ^ (Case III). (5.49) 

Similar constructions can be given in Cases I and II. For instance, we may build the first 
row of Na in terms of a meromorphic differential uj with prescribed simple poles and 
residues as shown in the following tables: 

(5.50) 



X 


P 


9 0i,2 


03,4 


001,2 


003,4 


Res(a;, x) 


-1/2 


-1/2 -1/2 


-1/2 


1/2 


3/2 



X 


P 


q 0i,2 




001,2 


003,4 


O3 O4 


Res(a;, x) 


-1/2 


-1/2 -1/2 


-1/2 


1/2 


3/2 


-a/2 a/2 



for Case I, and 

Ti n Hi n — ri nm n nm a On Oa 

(5.51) 

for Case II. Such constructions will again give us ()5.48p . but the behavior near the origin 
is different: 

Naiz) = 0(z-^/^), as z ^ (Case I), (5.52) 
7V„(z)diag(z^/^z^/^z"/^z-"/2) = 0(1)^ asz^O (Case II). (5.53) 

In what follows, we will construct local parametrices near each branch point. 
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5.5 Parametrix near the nonzero branch points 

Near each of the branch points in {p, — ri, — ra} \ {0}, we need to build a local 
parametrix using the 2x2 RH problem for Airy functions. This construction is 

very well-known and we omit the details. 



5.6 Parametrix near the branch point (hard edge) 

The local parametrix at the origin is built by means of modified Bessel functions. 
Such a construction was given before for a 2 x 2 RH problem in [39] and in a 3 x 3 setting 
in [37, 42j. Our construction will be similar to the one in ^37j. For — 1 < a < there are 
some extra complications that were solved in an ad hoc way in [37]. We will be facing 
similar problems and solve them in a more conceptual way. 



5.6.1 RH problem for modified Bessel functions 

First we recall the construction of [35]. Define the 2x2 matrix valued function 
^Bessel (-^^ by 



^Bessel ^^^j 

for I arg C| < 27r/3, 



/«(2CV2) ^/^„(2Ci/2) \ 



^Besse... ^^^^ (2(-C)^/^) (2(-C)^/^) 

for 27r/3 < argC < vr, and 



{H^a^) (2(-C)^/') {H^a'>) (2(-C)^/'), 



(i) 
a -I 



for — vr < argi^ < — 27r/3, where la and Ka are modified Bessel functions and H, 
z = 1, 2 is the Hankel function; cf. [U Chapter 9]. 

Denote with jj, j = 1, 2, 3 the complex rays {( £ C \ arg C = (j + l)'''"/^}. According 
to [39], satisfies the following RH problem: 

(1) ^Bessel analytic in C \ Uj=i 7j- 
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(2) With the rays 7j, j = 1, 2, 3 ah oriented towards the origin, jumps 



^Bessel ^Bessel ^ ^ 




on 71, 



on 72, 



on 73. 



(3) Uniformly for — t- 00 we have 



(5.54) 



(4) As C — in I argCI < 27r/3 we have 



Bessel 



(C) 




, if a > 0, 

if a = 0, 
if Q < 0. 



(5.55) 



5.6.2 Construction for a > 

In the construction of the local parametrix at the origin we will restrict ourselves to 
Case II, which is in some sense the most general and difficult of the three cases. The 
modifications for Cases I and III will be briefly commented at the end of this section. 

We denote by Bs a small fixed disk with radius 5 > centered at the origin, such 
that it does not contain any other branch points and consider all the jump matrices Jt 
of T restricted in Bg; cf. item (2) in Proposition 15.41 We note that the (2,3) entry of 
Jt on (0, 6) is e-2"<^2,+ . gince Re (/>2,+ = Re (A2 - A3) > c> on (0, 6) (cf. LemmaiZD, 
this entry is exponentially small as n — )• 00. Similarly the (4, 3) entry of Jt on (—J, 0) is 
_^2n<f)3 oigQ exponentially small. 

By neglecting the exponentially small entries in the jump matrices for T{z) in Bs, 
we are led to the following RH problem for a 4 x 4 matrix valued function Q{z): 



(1) Q is analytic in Bs \ (Ai U 
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Figure 7: Contours for the local parametrix around in Case II, for a > (left picture) 
and for —1 < a < (right picture). 

(2) On each of the oriented contours in the left picture of Figure [71 Q has a jump 
Q+ = Q-Jq, where 



Jq = ' 



diag 11° j ,e-^™,e*™ | , on Ai n 5^ = (-6,0), 



(3) As z — )• outside the lens enclosed by A^ we have 

Q{z) diag(|z|"/2, |z|-°/2) = 0(1), if a > 0, 

Q{z) diag((log \z\r\l, (log \z\)-\ 1) = 0(1), if a = 0. 



(5.56) 



(4) As n — > 00, 



Q(z)=iV„(z)(/ + 0(-) ), (5.57) 



n 



uniformly for z £ dBs \ (Ai U A^), where Na is the global parametrix built in 
Section 15. 4[ 

Before solving this RH problem, let us consider the function 

h{z) := {M^)f. (5.58) 

Since the density of pi — fxi blows up like 0(1x1^^^"^) as x — )• 0, we obtain from (j4.14p 
and the definition of (piiz) in ()4.38p that fi{z) is analytic in Bg and gives a conformal 
mapping from a neighborhood of the origin onto itself, such that fi{x) is real and positive 
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for X G (0,(5). For our purpose, we may deform the contours A]*^ near such that /i 
maps n Bs to the rays with angles 27r/3 and — 27r/3, respectively. 
We now define 

Q{z) = E{z) diag (ai*^^^^^'(n2/i(^))cTi, diag(e"<^i(^), -e-"<^^(^), 1, 1), 

(5.59) 

/o l\ 

where o"i = I I , and where the prefactor E{z) is analytic in and is chosen to 

satisfy the matching condition on dBs, see below. With this definition, the items (1), 
(2) and (3) in the RH problem for Q{z) are satisfied, by virtue of items (1), (2) and (4) 
in the RH problem for \I/Bessei^^-) 

To achieve the matching condition in item (4) of the RH problem for Q{z), we take 
E{z) in (i539]l as 

i^(z) = iV,(z)diag(l,-l,l,l)diag 4^ I ' 




V2 

X diag((27rn)-i/2/i(z)-V4, (2^n) V2/l(z)l/^ 1, 1). (5.60) 

Obviously E{z) is analytic for z £ Bs\ K^. Moreover, one checks that E{z) is also 
analytic across (—5,0). Finally, since /i(z)^/'^ = 0{z^^^) as z — )■ and using ()5.49p we 
find that 

^(z) = 0(z-i/2)^ ^^0, (5.61) 

so E(z) cannot have a pole at zero. We conclude that E(z) is analytic in the disk Bg. 
By virtue of (|5.54p . the matching condition ()5.57p in the RH problem for Q is satisfied. 
In summary, we have established the following proposition. 

Proposition 5.5. The matrix valued function Q defined by (|5.59p and (|5.60p satisfies 
the conditions (l)-(4) of the RH problem for Q . 

In view of the fact that det Na{z) = 1 and det xjr Vessel _ ^ ^g^g [3S]), we see that 

detQ(z) = l, zeBs. (5.62) 
If we take Q as the local parametrix for T, the final transformation would be defined as 

R{z) = T{z)Q{z)~^, z G Bs. (5.63) 
Then R would be analytic in Bg \ (0, oo) with the following jump on (0, 5), 
R-{x)-^R+{x) = Q{x)T^{x)-^T+{x)Q{x)-^ 

= Q{x){I + e-2"'^2'+(")S23)Q(2:)-' 
= / + e-2"*2,+(-)Q(3;)^23Q(x)-^ (5.64) 
Similarly, the jump on (—5, 0) would be 

R-{x)-^R+{x) = 1- e-*™e2"'^^(^')Q+(x)-E43Q+(x)~^ (5.65) 
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Lemma 5.6. For a > 0, the matrix Q{x)E23,Q{x) is bounded as x — t- 0, x > 0. 

Similarly, Q^{x)E43Q-^-{x)^^ is bounded as x — )• 0, x < 0. 

Proof. We prove the first statement only; the second fohows similarly. Observe that 

Qix)E23Qix)-^ = Qix) (o 1 o) fo 1 o)Q{x)-\ (5.66) 



If Q > 0, it follows from ()5.56p and det Q{z) = 1 that both Q(x) (o 1 o) and 



1 oj Q{x) ^ behave like 0(z"/^), and the lemma follows. 
If a = 0, we use (l535]l and det ^'Bessei(^) ^ ^ 

^B„A-'(^) / OW 0(1) \ ^^^^ 
' \o(log|C|) 0(log|C|)/ 

Inserting this formula in (|5.59p . we have 



-1/ 



[X 




diag(( ^^^^ \ ,0{1),0{1)\ E-\x), (5.67) 

O(logx) 0(logx 



as X —7- and x > 0, where £'~^(x) is bounded near the origin. Since Q(x)(0 1 0)^ 
is bounded by (|5.56p and (0 1 0)Q{x)^^ is bounded by ()5.67p . the lemma follows for 
a = as well. □ 

The above lemma, together with (|5.64p and the fact that Re (j)2,+{x) > c > for 
some constant c on (0, 6), implies that, if a > 0, the jump matrix for R is exponentially 
close to the identity matrix as n — )• oo, uniformly for x G (0,(5). A similar statement 
holds for the jump on {—5, 0). We therefore take the parametrix P = Q in case a > 0. 

5.6.3 Construction for — 1 < a < 

If — 1 < a < 0, then the (2,3) entry of the jump matrix Jt on (0,5) and the (4,3) 
entry of Jt on {—5, 0) cannot be simplify neglected. Hence, we need to build a 4 x 4 
matrix valued function P such that 

(1) P is analytic in Bs\{RU Af). 

(2) For each of the oriented contours shown in the right picture of Figure [71 P has a 
jump P+ = P-Jp, where 



Jp 








on Ai n ^5 = (-(5,0), 

on Af n Bs, 
on (0,(5). 
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(3) P{z) behaves near the origin hke: 

P(z) = 0(z-/2), 

(4) As n 



0{z 



a/2^ 



OO, 



P(z) = iV„(z) / + 



n 



(5.68) 



uniformly for z £ dBs \ (Ai U A^), where Na is the global parametrix built in 
Section [5l 



We use the matrix valued function Q given by ()5.59p and ()5.60p . which works as the 
parametrix for the case a > 0, to construct P. More precisely, recall the sectors Oj, 
J = 1, . . . , 4, in the right picture of Figure [3 By setting 

S{z) = A-^{z)Sj{z)A{z), 



z £ Qj, 



(5.69) 



where A is defined in ()5.2ip and Sj, j = 1, ... ,4 is a constant matrix that depends on 



the sector fij, we construct P in the form 



Piz) = Qiz)Siz). 



(5.70) 



The construction of these four matrices Sj is different in the three Cases I-III; we will 
again explain it for Case II. 

Let us give names to the jump matrix Jp on the four jump contours around zero as 
follows. Jp^i is the jump matrix of P on (0, 5), Jp^2 is the jump matrix on A^, Jp^3 on 
(—(5, 0) and Jp_4 on A^,". Similarly we define the jump matrices Jg^fc, /c = 1, . . . , 4. Note 
that Jq_2 = Jp,2, Jq,a = Jp,A, Jq,i differs from Jp^i in the (2,3) entry, and Jq_3 differs 
from Jp^3 in the (4, 3) entry, 

To obtain the correct jumps for P, the matrices Sj in ()5.70p should be such that 

AJp}A-i = 5f ^(AJ-\A-i)54, 
AJp,fcA-i = 5^-i(AjQ,fcA-i)5fc_i, for A: = 2,3,4. (5.71) 

Multiplying these equations into a telescoping product form yields the following relation 
for 5i: 

A{Jp\jpAJp,3Jp,2)A-' = S^'A (Jq\JqaJq,3Jq,2) A~'Si. (5.72) 



Now A( Jp^ Jp^4 Jp^3 Jp^2)A ""^ and A( Jq"*]^ Jq^4 Jq^3 Jq^2)A ^ are both constant matrices. A 
calculation shows that they have the spectral decompositions 



A{Jp\jp,AJp,'iJp,2)A 



-1 







-1 



V 

diag(e' 








\ 



/ 

- E' 



VpDVp\ (5.73) 



21 



VqDVq 



(5.74) 
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with common eigenvalue matrix 



D = diag(e-*™,e'™,e-'''",e*''"), (5.75) 
and where the eigenvector matrices Vp and Vq can be chosen as 

Vp = I+ / {E21 + e-'™ii;23 + i?43), Vq = I+ / E21. 
z^sln(7^Q^) 22Sin(7raj 

On account of the spectral decompositions (|5.73|) - (j5.74|) . the equation (|5.72|) can be 
rewritten as D{Vq'^ SiVp) = (Vq ^SiVp)/?, which allows the solution Si = VqVp^. So 
we can take 

Si = VqVp^ = I- ^. A e-'^^E2^ + i?43). (5.76) 

2zsm(7raj 

The constant matrices 5*2, 5*3 and ^4 are then obtained recursively from Si and ()5.7ip . 
which now read 

52 = I + 77— (^13 - e-'^^E2^ - i^43) , (5.77) 

53 = /-T7— V^(^13 + e*™^23 + ^43), (5.78) 

54 = / - -— (e^™i?23 + i^43). (5.79) 

2«sm(7ra) 

With S'fe, /c = 1, . . . ,4, given in (j5:76|) - ([579|) . we conclude that P defined in (fSTTO]) 
indeed satisfies the items (l)-(4) in the RH problem for P{z) stated above. 

The above description was for Case II. For Cases I and III similar constructions can 
be given. The key point is that in each case, the cyclic products of jump matrices in 
the left hand sides of (j5.73p and (j5.74p . respectively, allow spectral decompositions with 
common eigenvalue matrix (|5.75p . We leave these constructions to interested readers. 

Combining the results in this subsection and Section 15.6.21 we take P = QS as 
the parametrix in for Case II, where 5" = / if a > 0, and S is given by (|5.69p . if 
-1< a < 0. 

5.7 Final transformation 

Using the local parametrices P, P^^^y and the global parametrix N^, we define the 
final transformation as follows 



R{z) 



T{z){P^"'y (z)) ^, in the disks around {p, q, — ri, —r^} \ {0}, 
T{z){P(z))^^ , in the disk Bs around the origin, (5.80) 
T{z){Na{z))^^ , elsewhere. 



From our construction of the parametrices, it follows that R satisfies the following 
RH problem: 
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Figure 8: Jump contours for the RH problem for R in Case I, for a > (top) and 
for —1 < a < (bottom). 



(1) i? is analytic in C \ E^, where S/j depends on a and is different for all three cases. 
For Cases I and III, an illustration is shown in Figures [8] and [9l respectively. 



(2) R has jumps R+ = R-Jr on S/j that satisfy 

jR{z)=I + 0{l/n), 
uniformly for z on the boundaries of the disks; 

J,j(x) = / + 0(x"e-^"), 



(5.81) 



(5.82) 



on (0,(5) (Case I), or on (—5,0) U (0,(5) (Case II), or on (—(5,0) (Case III), for some 
constant c > 0, if (a; > 0; 

JH(z) = / + 0(e-^"l^l), (5.83) 
on the other parts of S/j, for some constant c > 0. 

(3) R{z) = 1 + 0{l/z) as z ^ oo. 

Then, as in [151 [13 EZ]) we may conclude that 



Riz) = I + 



1 



n{\z\ + 1), 

as n — 7- oo, uniformly for z in the complex plane outside of S/j. 



(5.84) 



6 Proof of Theorem 12.3 



With the above preparations, we are ready to establish Theorem 12.31 For this pur- 
pose, we need to represent the correlation kernel Kn{x, y) in terms of T by following the 
sequence of transformations Y ^ X ^ T m. the RH analysis. 
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Figure 9: Jump contours for the RH problem for R in Case III, for a > (top) and 
for —1 < a < (bottom). 



Let x,y £ A2. Recalling the representation of Kn{x,y) in (|2.34|) . we first find from 
(15:711 that 



Knix,y) 



2'Ki{x — y) 
In view of ()5.18p . this becomes 



Kn{x,y) = \ (0 e"^3,-(?/) U^\y)U+ix) (0 e'^^^.+W q o)' • 



e"— 



2iTi{x — y) 
By (lOHjl . it follows that 



Kn{x,y) 



27^i{x-y)\^ -e-^^.-(^) e"^3,-(.) 0) r-Hy)r+(x) 



Now, we obtain from ()5.84p and standard arguments (e.g. [6, Section 9]) that 

T^\y)T+{x) = I + 0{x-y), as x ^ y, (6.2) 



uniformly in n. Also note that Xs^± = A2,=p. Thus, by taking y — )• x, it follows from ()6.2p 
and L'Hopital's rule that 

Knix, x) = ^(e2,+(x) - e2,-(x)) + 0(1) = n^(x) + 0(1), (6.3) 
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where the last equahty follows from (|4.15|) . and so 

lim —Kn(x,x) = -T^(x), X £ A2. (6.4) 

n-s>oo n ax 

If X € M"*" \ A2, a similar argument as given above yields 

lim -Kn{x,x) = 0. (6.5) 

n— !-oo n 

This completes the proof of Theorem 12.31 
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